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Abstract. We discuss the physical consequences of a duality between two models 
with quenched disorder, in which particles propagate in one dimension among random 

■ traps or across random barriers. We derive an exact relation between their diffusion 
C/3 , fronts at fixed disorder, and deduce from this that their disorder-averaged diffusion 

fronts are exactly equal. We use effective dynamics schemes to isolate the different 
physical processes by which particles propagate in the models and discuss how the 
duality arises from a correspondence between the rates for these different processes. 



1. Introduction 



(N 

I> How do classical particles in heterogeneous environments propagate? This question 

arises in many contexts, from glass- forming liquids and colloids [H |2] , to biomolecules 
^ I moving in the crowded environment of the cell [3j, to electrical properties of disordered 
materials [1] . Diffusion is the most familiar mode of propagation, but if the heterogeneity 
involves a sufficiently broad distribution of time scales then slower motion (subdiffusion) 
is possible: for some reviews of this well-studied field, see [5], El 13 |8] . 

Important experimental evidence of subdiffusive motion in disordered environments 
was obtained a long time ago, through the electrical conductivity of hoUandite (a quasi 



' one-dimensional ionic conductor). Here, the charge-carriers can be modelled by classical 
particles on one-dimensional chains: their anomalous low frequency conductance was 
attributed to subdiffusive motion, and modelled by a simple one-dimensional disordered 
model [Ij. More recently, subdiffusive propagation has been important when considering 
dynamically heterogeneous relaxation in glass-formers. There is considerable evidence 
(for example, [H [21 [9l [10]) that while time scales grow very fast with decreasing 
temperature, the length scales associated with dynamical heterogeneity increase rather 
slowly. In the picture of the glass transition based on dynamical facilitation pTl[T2l lT3]. 
this implies a subdiffusive propagation of mobility through the system [H]. Several 
possible mechanisms have been proposed for such an effect, such as directional kinetic 
constraints JT2\, or a strong dependence of the dynamics on the local free- volume |17j . 
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Thus, understanding how subdiffusive motion arises in simple models remains an 
important challenge. 

In this article, we discuss the physical processes that lead to subdiffusive 
propagation in two very simple one- dimensional models. They are directly relevant 
for materials such as hollandite [1], and for the models of glassy behaviour discussed 
in [17], but they have also been related to the motion of defects in disordered magnets, 
to disordered elastic chains and to networks of resistors and capacitors (see O [7] for 
reviews). In the first model, diffusive motion is frustrated by the presence of sites where 
the particle gets trapped for a long time; in the second, there are barriers between sites 
across which motion is very slow. These models are related by a duality relation between 
their master operators. This relation has been noted before [5[ [71 [18] , but with limited 
discussion of its physical content. Here, we revisit the duality between the two models, 
focussing on their diffusion fronts, and derive exact relations between these which do not 
appear to have been noticed previously. We connect the evolution of the diffusion fronts 
with the events in which the particle escapes from deep traps or surmounts large barriers. 
We elucidate the physical processes underlying the duality, and use it to motivate an 
effective dynamics for the model, following le Doussal, Monthus and Fisher p[9] . 

We define our models in Section [2] and discuss the duality relation in Section [31 We 
cast the duality relation for fixed disorder as a local relation between the propagators 
of our models, and show that it leads to equal disorder-averaged propagators in the two 
systems. In Section [H we describe the physical processes by which the models evolve, 
and use them to define effective dynamics schemes. The extent to which this dynamics 
captures the real-time evolution of the diffusion front is also discussed. We summarise 
our conclusions in Section [H 

2. Models 

We first define the random trap model [7]. Consider a single particle hopping on a 
chain of sites. A particle on site n may hop either to the left or right, and hops in both 
directions happen with the same (site-dependent) rate, Wn- The master equation for 
the distribution of particle positions, p^, is then 

{d/dt)pl{t) = Wn+ipl,^,{t) + Wn-ipl.,{t) - 2Wnpl{t) (l) 

The superscript 'T' distinguishes quantities for the trap model from those for the barrier 
model. In the latter [5l [7], the rates Wn are associated not with sites, but with links 
between them. That is, hops from site n to n + 1 happen with rate Wn, as do hops from 
site n -|- 1 to n. In this case, the master equation is 

{d/dt)p^Xt) = wM^,{t) - p^{t)] + w„-i[pti(t) -p^it)] (2) 

In [7], it was noted that if the same set of rates Wn are used in both models, 
then the current in the barrier model = Wn{Pn — Pn+i) obeys the same differential 
equation as the rescaled probability density WnPn in the trap model. This duality 
between trap and barrier models was noted a long time ago by Dyson pSj: we refer 
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to it as 'trap-barrier duality'. As explained below, it follows that the eigenspectra of 
the two master operators are equal for these two models: this property is relevant for 
the impedances of linear chains of resistors and capacitors, and for the elasticity of 
random one- dimensional chains. However, in understanding how particles propagate in 
disordered media, the most natural quantity is the propagator, or diffusion front, which 
we discuss in some detail below. 

It was further argued in [7J that trap-barrier duality can be understood by 
identifying the regions between large barriers in one model with the deep traps in the 
other. While very appealing intuitively, this argument is possibly somewhat misleading, 
in two opposite ways. On the one hand, it suggests that the duahty holds only after 
suitable coarse-graining to large lengths and times; we will see, however, that for the 
diffusion fronts it is valid on all length and time scales. On the other hand, viewing a 
barrier model as consisting of effective traps would suggest that, at least on large enough 
scales, the models will behave effectively identically. This also is not quite correct. For 
example, the barrier model has a stationary distribution in which the particle occupies 
each site with equal probability, while the trap model exhibits aging behaviour if it is 
initialised in a uniform state [20]. Quantitatively, consider a particle that starts from 
the origin at time zero, and measure its mean square displacement between times t' and 
t > t'. Averaging over the random rates Wn according to a translationally invariant 
distribution, we have (see Section [3]) that 

R\t,t') = {[xit)-xit'W) (3) 

depends only on the time difference t — t' in the barrier model (the angle brackets 
represent an average over the stochastic dynamics, and the overbar represents the 
average over the rates Wn)- On the other hand, in the trap model, R^{t,t') generically 
depends on both time arguments, and for the disorder distributions we will consider it 
never reaches a stationary regime in the thermodynamic limit. 

Bearing this and other differences in mind, we now use the relationship between ([T]) 
and ([2]) to investigate the extent to which the barrier model can be modelled by a set 
of effective traps. We establish relationships between the propagators (diffusion fronts) 
of the two models, and use this analysis to develop a physical picture of propagation 
in these systems. We also show that the duality implies not just that the long-time 
behaviour of these models are the same, but that their disorder averaged propagators 
are identical at all times. We contrast the extent to which the models are similar with 
differences between them, including their aging behaviour and the properties of the 
typical particle trajectories. 

3. Duality 

It is convenient to define the master operator for the barrier model, , in terms of 
its matrix elements W^^, by writing ([2]) as {d/dt)p^{t) = —J^m^wmPmi't)- We specify 
periodic boundaries on a chain of length L, so this operator has L eigenvectors q^. The 



Duality between random trap and barrier models 



4 



eigenvector with zero eigenvalue gives the steady state q^''^'^ = 1/L. The dynamics of the 
model obeys detailed balance with respect to this trivial distribution, as is clear from 
the symmetry = W^^. For the trap model we call the master operator W^; it too 
obeys detailed balance but with respect to the non- uniform steady state q'^''^'^ oc W~^. 

The expression of trap-barrier duality in terms of propagators can be derived from 
the simple fact that, if Pn{t) is a solution of the barrier model master equation ([2]) then 
Pn{t) — Pn+iit) is a solution of the trap model equation ([1]). Applying this to a barrier 
model eigensolution p^it) = exp(— At), it follows that = q^ — qn+i is an eigenvector 
for the trap model master operator with the same eigenvalue A. This is a one-to-one 
relation between the eigenvectors with nonzero eigenvalues of the two master operators, 
and so their spectra are identical as claimed in Section [2l (Inverting the differencing 
relation to get from q"^ back to q^ in principle gives an undetermined constant, but this 
is fixed by the requirement that nonzero eigenvectors obey J2n1n = 0-) The exception 
is the steady state, where differencing the barrier model steady state q^'^"^ = 1/L gives 
a vanishing result rather than the trap model steady state. 

The propagator or diffusion front of the barrier model, G^^^(t), is defined as the 
solution of ([2]) with initial condition pJ^(O) = 6n,i- This can also be obtained from matrix 
elements of the time evolution operator: G^^^it) = [exp{—W^t)]ni. The symmetry of 
the barrier master operator then implies that also the propagator is symmetric 
under interchange of arrival and departure sites 

G^nUt) = Gtnit) (4) 

as required for a system obeying detailed balance with respect to a uniform distribution. 

We now turn to the trap-barrier duality relation for the propagators. Since 
Pnit) = G^^j(t) is the solution to ([2]) with initial condition pJ^(O) = Sn,i, it follows 
that Pnit) = G^^^{t) — G'^+i^.j(t) is a solution of ([1]), with initial condition p^{0) = 
Sn,i — Sn+i,i = Sn,i — ^n.i-i- Defining the propagator for the trap model G'^^j(t) as the 
solution to ([1]) with initial condition p^(0) = 6n,i, it follows (by linearity of the master 
equation) that 

Gn^iit) - G'^+i^i(t) = Gj^^i{t) - Gj^^i_i{t) (5) 

This identity, which holds for all realisations of the disorder {W^n}, is our desired exact 
statement of the trap-barrier duality relation, in terms of the propagators of the two 
models. It can be used to express the trap model propagator in terms of the barrier 
model one, and vice versa. For example, since the diffusion front vanishes at large 
distances, we have 

GfUt) = EiGfU^) - Gf^.Ut)] = EiGjUt) - GU-M- (6) 

j=n j=n 

(We require that the chain length L is large enough that G^^-j does indeed vanish at 
large j. That is, L should be taken to infinity before taking any limit of large time.) 
A corresponding relation can be used to express G^ in terms of G^. We show typical 
propagators in Fig. [1], and give a geometrical interpretation of ([5]): the difference of the 
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Figure 1. (Top) Sample propagators, Gj^^i and Gj^^j, as a function of final site n, 
for trap and barrier models at /i = 3. The set of disordered rates is the same in each 
case, and we take i = and t = 2^^. Small rates lead to sites with peaks in the trap 
model propagator, and to large changes between adjacent sites in the barrier model. 
(Middle) Propagator maps for the same systems: the propagators at fixed time are 
plotted as functions of both initial and final positions. The dashed lines indicate n = i 
and the dotted lines with arrows indicate the slices for which we plot data in the other 
panels of this figure. In the trap model, the probability is concentrated on final sites 
with small escape rates. In the barrier model, it is delocalised across 'effective traps', 
separated by barriers with small crossing rates. The duality relation ([5]) states that 
the gradient of the trap model propagator with respect to i is equal to the negative 
gradient of the barrier model propagator with respect to n. The relevant gradient 
in the barrier model is localised on the upper and lower edges of the effective traps, 
and these are the final sites on which the probability is localised in the trap model. 
(Bottom) We plot the trap model propagator as a function of initial site i, for a fixed 
final site n = 4. This illustrates its gradient with respect to i. 
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barrier model propagator taken between neighbouring arrival sites equals the negative 
difference of the trap model propagator between neighbouring departure sites. 

In making the plots in Fig. [1] we considered the most common case where the Wn 
are independently and identically distributed (i.i.d.), with a distribution 

P{W) oc W^^/^^-\ 0<W <1 (7) 

With this choice, the long time behaviour of a particle is diffusive for /i < 1 and 
sub diffusive for fi > 1. More precisely, if a particle begins on a randomly chosen site, 
its displacement at large times scales as [TJ El] 

z = max(2, 1 + /i) (8) 

with logarithmic corrections for /i = 1. 

It is sometimes useful to interpret the rates Wn as arising from activation energies 
En associated with crossing a barrier or leaving a trap. We then write Wn = 
exp{—En/T), where T is the temperature and we have set Boltzmann's constant equal 
to unity. In this representation, assuming the distribution of the activation energies is 
P{E) = exp{—E) with E > 0, we can identify fi with 1/T. 

We now consider averaging the propagators and over realisations of the 
disorder. For any distribution of the disorder which is translationally invariant (which 
includes all choices in which the Wn are i.i.d.), the disorder-averaged propagators depend 
only on the difference m = n — i. Using a bar to denote the disorder average as before, 
we thus have from the duahty relation (l5l) 

Gm{t) - = G^{t) - G'„+i(t). (9) 

Since both propagators must vanish as |m| ^ oo, it follows that 

Gl{t)=Gl{t) (10) 

This exact coincidence of the disorder averaged diffusion fronts is a surprising result. 
For a given realisation of the disorder, the propagators for the two models are very 
different, as shown in Fig. [TJ For other quantities the difference is even more striking. 
Consider for example the quantity R'^{t,t') defined in ([3]) above. Since all propagators 
decay to zero at large distances, we can replace the disorder average by an average over 
initial position i, so that R^{t,t') = {1/L) J2k,n{^n — Xk^Gn^kit — t') J2iGk^i{t'), where 
Xn is the position of site n. But now the symmetry (j4]) of the barrier model propagator 
together with conservation of probability implies Y.iGk^i{t') = J2iGi^k{t') = 1- Thus 
R^{t,t') depends on t and t' only through the propagator Gn^mit — t'), and therefore 
only through the time difference t — t'. The argument evidently generalizes to the entire 
(disorder-averaged) distribution of x{t) — x{t'). This time translation invariance arises 
because a particle moving among barriers is equally likely to be on any site, whether 
or not it is bordered by a high barrier. In the trap model, on the other hand, the 
particle spends (for /i > 1) most of its time in deep traps whose escape rate decreases 
with the age t' as ~ [7J. To arrive at the asymptotic scaling of R'^(t,t') 

in the trap model, we take t > t', and note that moves to left and right are always 
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equally likely, so that {[x(t) — x(t')]x(t')) = (for any realisation of the disorder, as 
long as the system size L is taken to infinity before any limit of large time [22]). Thus, 
R'^{t,t') = {x^{t)) — {x'^{t')), which for large t and t' scales as t^/^ — t'^/^. Clearly, time 
translational invariance is obtained only in the diffusive case {z = 2 ot fi < 1). 

The difference between the two models in the subdiffusive regime /i > 1 can also 
be seen at the level of typical trajectories. For example, let N be the mean number of 
hops made by a particle up to time t. Time translation invariance in the barrier model 
implies that ~ t. To get the scaling ([H]), the displacement must then grow with 
the number of hops as \x\ ~ A^^/^, i.e. more slowly than for simple diffusion. On the 
other hand, in the trap model, the typical number of hops associated with trajectories 
of displacement scales as in simple diffusion: ~ because escapes from any 
trap occur with equal probability to the left and the right. This then implies a number 
of hops growing only sublinearly in time, N ~ t^/^ [7]. 

To summarize thus far, the trap and barrier models have on average equal 
diffusion fronts, but many other physical properties differ. Motivated by this surprising 
observation, we now discuss the duality relation at fixed disorder in more detail, in order 
to understand how the equality of the disorder-averaged diffusion fronts arises. 

4. Effective dynamics 

In Fig. [1] we showed the propagators for trap and barrier models, for a single realisation 
of the disorder (with /i = 3) and at a given time. To understand how the propagators 
evolve in time we exploit an 'effective dynamics picture' in the spirit of [19]. A very 
similar scheme was applied to the trap model in [23l [2^ . The effective dynamics is 
based on an assumption of well-separated time scales, which is valid in the limit of large 
yU taken at fixed energies En- We first present the scheme for the barrier model. 

4.1. Barrier model 

To define the effective dynamics in this model, we separate those barriers which are 
relevant on a time scale t from those which are irrelevant. Large barriers (those with 
small transmission rates W) tend to be relevant, because they limit the motion of the 
particle. Once we have identified a set of relevant barriers, we assume that particles 
move rapidly between them, but never cross them. Thus, for an initial site m, we have 

- 1 0, otherwise 

where i is the index of the nearest relevant barrier to the left of site m, and j the 
index of the nearest relevant barrier to its right (see Fig. [2]). These two barriers 
define the 'effective trap' within which the particle is localised. In the propagator map 
representation of Fig. [1], the effective dynamics models the propagator as a series of 
non-overlapping blocks. As the set of relevant barriers evolves in time, the map evolves 
by instantaneous events in which two adjacent blocks coalesce into a single larger one. 
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Figure 2. (Left) Sketch illustrating a single step of the effective dynamics for the 
barrier model. Barriers i, j and k are all relevant at time to, so G'J^^^(to) is localised 
between barriers i and j. At the later time ti, barrier j has become irrelevant 
{Kjti > 1). (Right) A sketch of the corresponding behaviour for the trap model. 
The propagator is localised on sites i and j at the initial time Iq and on sites i and k 
at the later time ti. The value of the propagator on these sites is in inverse proportion 
to their distance from the initial site m. 



To obtain the time scale on which barriers cease to be relevant, suppose that the two 
barriers delimiting the effective trap have well-separated transmission rates, Wj ^ Wi. 
On time scales much smaller than the diffusion front evolves in time only by 

transmission through barrier j. The diffusion front increases within a neighbouring 
effective trap, which is formed by barrier j, and the nearest relevant barrier to its right 
(let the index of this barrier be k, as in Fig. [2]). We now further assume that the motion 
between barriers j and k is fast compared to transmission through them. In that case, 
we can approximate the propagator by 

i Pi{t), J <n<k 

where po and pi are probability densities within the two effective traps (to the left and 
right of barrier j: see Fig. |2]). The probabilities in these traps are po/o, and pi/i, where 
Iq = j — i and li = k — j are the sizes of the effective traps. These probabilities evolve 
as dt[lopo{t)] = -dt[lipi{t)] = Wj[pi{t) - po{t)]. Hence, 

dt[po{t) - pi(t)] = -K,[po{t) - p^{t)], (13) 

and we identify 

as the rate for equilibration between these two effective traps. 

Still working within our assumptions of well-separated time scales, we can use this 
procedure to define the time evolution of the set of relevant barriers. At small times, all 
barriers are relevant. For a given set of relevant barriers, we characterise barrier j by its 
value of Kj, where the lengths Iq and li are the distances to the nearest relevant barriers 
to the left and right of barrier j. (When a relevant barrier is removed, a new effective trap 
is formed, and the values of Kj on the adjacent barriers change according to the width 
of the new trap.) Thus the effective dynamics evolves in time by successive removal of 
the barrier with the largest value of Kj. To obtain the set of relevant barriers at time 
t, this procedure is iterated until all relevant barriers have Kjt < 1. The propagator is 
then constructed using (fTTl) . as described above. 
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Figure 3. (Left) Propagator map for the barrier model at /i = 10 and t = 2^^. 
This was calculated by histogramming the corresponding trap model propagator 
and using (Right) Propagator map for the effective dynamics, using the same 

realisation of the disorder and the same final time. The agreement is quite good, 
although the real propagator has some fine structure not captured by the effective 
dynamics. 



This approach is similar to that of [T9], differing only in that the sizes of the 
effective traps enter the transition rates. This is necessary for this model, since these 
factors affect the scaling of the diffusion front (see Section 1^73]) . We also note that the 
rate Kj depends symmetrically on the width of the effective traps to the left and the 
right of the barrier. This ensures that the effective dynamics for the diffusion front 
preserves its symmetry (jl]). 

We now revisit our assumptions regarding well-separated time scales. These 
assumptions were (i) that the two barriers delimiting each effective trap have well- 
separated rates, allowing us to treat transmission over them independently, (ii) the 
rates associated with successive stages of the effective dynamics are well-separated, and 
so the time taken for a given barrier to become irrelevant is well-approximated by K~'^, 
and (iii) equilibration within each of two neighbouring effective traps is faster than 
transitions between them [so that (fT2|) can be applied]. Writing Wn = exp{—En/T) as 
before, with T = 1/ fi, we recall that the distribution of energies is P{En) = exp{—En)- 
For the first assumption, the time scales for crossing adjacent barriers are well-separated 
if the difference between their two energies is much greater than T. In the limit of low 
temperature (or large /i) the probability of finding two adjacent barriers with rates 
within T of each other vanishes, so this assumption holds. For the second assumption, 
the same argument implies that the rate for any stage of the dynamics is well-separated 
from the (larger) rate of the previous stage. Finally, for the third assumption, the rate 
for equilibration within an effective trap is of the order of the rate of a previous stage 
of the effective dynamics, which is well-separated from the rate of the current stage as 
per the second assumption. 

In Fig. [3l we compare the propagator for a given realisation of the disorder with the 
prediction of the effective dynamics, at = 10. The agreement is quite good, justifying 
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our assumptions of well-separated time scales at this value of fi. As /i is decreased from 
this large value towards unity, these assumptions break down. The result is that the 
propagator cannot be well-approximated by non-overlapping blocks in the density map 
representation (recall the structure of the propagator at /i = 3, shown in Fig. [1]). We 
discuss the agreement between effective and real dynamics in more detail in Section 14.31 
at the level of the disorder- averaged diffusion front. First, though, we construct an 
effective dynamics for the trap model, related to the effective dynamics of the barrier 
model by the duahty relation ((51). 

4-2. Trap model 

Since we have the duality relation ([5]), the effective dynamics that we have described for 
the propagator of the barrier model has a corresponding effective dynamics in the trap 
model. We begin by describing this effective dynamics in terms of physical processes in 
the trap model, and then we show that it is indeed the dual of that of the barrier model. 

We assume that for a given initial site m, the propagator is concentrated on two 
'relevant' traps on sites i and j, with i < j. As before, if time scales are well-separated 
then we can consider these two sites separately. If Wi <C Wj, the relevant physical 
processes are those in which the particle escapes from trap j and falls into another 
relevant trap: either the one on site i, or one located on a site k > j (see Fig. [2]). 
The particle hops from site j to j ' + 1 with rate Wj. However, the probability of it 
reaching site k before being reabsorbed on site j is l/{k — j). (Particles always hop left 
and right with equal probability, so this combinatorial factor is the same as for simple 
diffusion.) Thus, the rates for motion from site j to sites i and k are Wj/{j — i) and 
Wj/{k — j) respectively. We therefore identify the total rate for moving from site j to 
either neighbouring trap as Kj, as defined in (IT^ . This is the time scale on which the 
trap on site j becomes irrelevant. Hence, the effective dynamics for the set of relevant 
rates in the trap model is the same as the effective dynamics for those in the barrier 
model. 

It remains to find the value of the propagator on these relevant sites. Starting from 
a single trap on site m, with i < m < j, the rates for arriving on sites i and j are in 
inverse proportion to their relative distances from m. Therefore, if at time to trap m has 
become irrelevant then G^^(to) = (j —tti)/ (j — i) and G'J<_„(to) = {jn — i)/ (j —i). In 
the next step of the process, site j becomes irrelevant, by time ti say, and the probability 
on that site is redistributed onto sites k and i, with probabilities inversely proportional 
to their distances from site j. Thus, we can relate the propagator at ti to that at to by 
G]^^{ti) = Gj^^{to) + ^Gj^rnito) = ^rf and Gj^j^{ti) = |5fGj^^(to) = fff • The 
relative probabilities of the two relevant sites at time ti are thus in inverse proportion to 
their distances from the initial site m, just as they were at time to. Hence, this property 
is maintained as successive barriers become irrelevant. 

These arguments motivate us to define our effective dynamics for the diffusion front 
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as: 




J — m 



m — I 



(15) 



I 



3 



where i and j are the nearest relevant traps to the left and right of the initial site m, 
and the set of relevant traps is determined by successively removing the trap with the 
smallest value of Kj. A similar result was derived in [23], although relevant traps were 
identified by their values of Wj in that work. These choices seem to be equivalent in 
the limit of large in which the effective dynamics schemes are valid. Returning to our 
analysis, it is then trivial to check that ( |T5l) and ( ITTl) are consistent with the duality 
relation ([5]), and that they are therefore dual to each other. The comparison between 
effective and real dynamics for the barrier model was shown in Fig. [3l Since both real 
and effective dynamics obey the duality relation, the agreement is similarly good for the 
trap model. 

Thus, we have identified the physical processes leading to motion in the trap model, 
and used these processes to motivate an effective dynamics for the diffusion front. This 
effective dynamics is related by the trap-barrier duality of ([5]) to that defined in the 
previous section for the barrier model. We trace this symmetry to the equality between 
the rates Kj for the (different) physical processes that occur in the two models. 

4-3. Disorder-averaged diffusion fronts 

Turning to disorder-averaged properties of the effective dynamics, the scaling ([8]) for 
/i > 1 arises directly from the form of the rate Kj given in (JT^ . To demonstrate this, we 
use the language of the barrier model. Let the typical distance between relevant barriers 
at time t be {i{t)). (Unlike distributions of times, all moments of the distribution 
of lengths are finite, so scahng arguments based on typical widths are valid.) Since 
relevant barriers have Kt < 1, the rates W associated with these barriers typically 
satisfy Wt < {(-{t)). For consistency, the typical distance between these barriers must 
scale as (i'(t)), so we have 



where P{W) was given in ([7]). This leads to t ~ £(t)^+^. All length scales in the effective 
dynamics scale with £(t), so this dynamics leads to z = 1 + /i, giving the correct scaling 
of the diffusion front for all /x > 1 [recall ([8])]. 

Following the reasoning in [19], scaling arguments can also be used to express the 
disorder-averaged diffusion front in terms of the fraction of effective traps of length ^ at 
time t. In the scenario studied in [191 EH], where the widths of all effective traps and 
the rates of the relevant barriers are all independent, this allows one to make significant 
progress. By contrast, in our effective dynamics the relevance of barriers is determined 
by Kj, and this introduces correlations between the rates and the widths of the effective 
traps. We have therefore not been able to determine the required distribution of trap 
lengths i analytically, and instead return to numerical simulations. 




(16) 
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Figure 4. (Top) Plots of {x'^{t)), comparing real dynamics (symbols) and effective 
dynamics (solid lines). The dashed lines show the power law scaling of ([5]), with 
unit prefactors (we plot (x^) = i^/(^+^)). The effective dynamics captures the time- 
dependence of the mean-squared displacement quantitatively, even at the relatively 
small value /i = 3. (Bottom) We show the disorder-averaged diffusion front G„(i) 
for real dynamics (points) and effective dynamics (lines) at fi — 10 (left) and /i = 3 
(right). The times are those of the final points in the top panel, t — 2**^ and t — 2^^: 
the data of that panel indicate that these times are long enough to be representative 
of the asymptotic large-time scaling. The effective dynamics agrees well with the real 
dynamics for /z = 10, but for the smaller /i the agreement deteriorates because the 
assumption of well-separated time scales breaks down. 



In Fig. m we compare the predictions of the effective dynamics with simulations 
of the trap and barrier models. At long times, we find good quantitative agreement at 
/i = 10. At yU = 3, the effective dynamics still captures the width of the distribution, 
although the absence of well-separated time scales leads to deviations in the tail. The 
diffusion front for the trap model was considered in [25], where it was predicted to 
decay as Gn{t) ~ exp[— for large n and /i close to unity (in our 

notation, the parameter /x of [25] is denoted by l/jj). While our results are consistent 
with this prediction, our statistics are not good enough to accurately determine the 
asymptotic scahng of the tails of the diffusion front. Overall, we would argue that our 
results indicate that the proposed effective dynamics does indeed capture the physical 
processes responsible for propagation in the trap and barrier models. 
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5. Conclusions 

We have discussed several implications of the duality between one-dimensional trap and 
barrier models, including ([5]), which is an exact local relation between their diffusion 
fronts valid for fixed disorder, and ( JTOl) . which states that their disorder-averaged 
propagators are equal. In addition, we have presented a unified effective dynamics 
for diffusion fronts in the large /i limit of the trap and barrier models. This dynamics 
captures the scaling of the diffusion front, and gives good quantitative estimates of its 
width. Our arguments indicate that it gives the exact shape of the diffusion front in 
the limit of extreme subdiffusion (/i oo), with agreement becoming less good as the 
system approaches the diffusive regime (/x < 1). 

The effective dynamics exploits the fact that in both models, the time scale on 
which traps and barriers become irrelevant depend on both their rates W, and on the 
spacings between them. In the barrier model, the dependence on spacings arises from 
the fact that if the relevant barriers are well-separated, the particle spends only a small 
fraction of its time adjacent to them. In the trap model, it originates from the fact 
that a particle leaving a given trap is likely to be reabsorbed there before it can reach a 
neighbouring deep trap. Despite the different physical mechanisms, these two processes 
affect the diffusion front in closely related ways, and mathematically result in the master 
operators for these processes having equal eigenspectra. 

Thus, for the purposes of the diffusion front, it is indeed appropriate to consider 
the barrier model as representing an environment of effective traps in which a particle 
diffuses. However, as we discussed, the trajectories by which the system evolves in each 
case are different. The effective dynamics illustrates the nature of these trajectories: 
in the barrier model, the particle diffuses rapidly within its effective trap, leading to a 
number of hops growing as ~ t. In the trap model, the particle is localised in deep 
traps, making excursions from its current trap which mostly lead to reabsorption within 
that trap. The mean hopping rate decreases with time as the system ages [20j, leading 
to ~ t^/^. It would be illuminating to investigate whether the duality relation ([5]) at 
fixed disorder can be cast as a relation between the particle trajectories in the barrier 
and trap models, to see whether these different scalings emerge naturally. 

Given the different physical processes underlying the dynamics of the models, 
we argue that the relationships between the diffusion fronts are, to a certain extent, 
coincidental. They do not represent a physical equivalence of the models themselves, and 
they do not generalise in a simple way to correlation functions other than the diffusion 
front. However, it seems to us that the equality of the disorder-averaged propagators 
and the associated algebraic structure might be exploited in further analytic work on 
these models. 
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